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Abstract: We present a detailed phenomenological study of direct photon production 
in association with a heavy-quark jet in pA collisions at the Relativistic Heavy Ion Col- 
lider (RHIC) and at the Large Hadron Collider (LHC) at next-to-leading order in QCD. 
The dominant contribution to the cross-section comes from the gluon-heavy-quark (gQ) 
initiated subprocess, making 7 + Q production a process very sensitive to both the gluon 
and the heavy-quark parton distribution functions (PDFs). Additionally, the RHIC and 
LHC experiments are probing complementary kinematic regions in the momentum fraction 
X2 carried by the target partons. Thus, the nuclear production ratio F(^^ can provide 
strong constraints, over a broad x-range, on the poorly determined nuclear parton distribu- 
tion functions which are extremely important for the interpretation of results in heavy-ion 
collisions. 
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1 Introduction 

Parton distribution functions (PDFs) are an essential component of any prediction involv- 
ing colliding hadrons. The PDFs are non-perturbative objects which have to be determined 
from experimental input and link theoretical perturbative QCD (pQCD) predictions to ob- 
servable phenomena at hadron colliders. In view of their importance, the proton PDFs 
have been a focus of long and dedicated global analyses performed by various groups; see 
e.g. Refs. [1-8] for some of the most recent studies. Over the last decade, global analyses 
of PDFs in nuclei - or nuclear PDFs (nPDFs) - have been performed by several groups: 
nCTEQ [9-11], nDS [12], EKS98 [13], EPS08/EPS09 [14, 15], and HKM/HKN [16-18] (for 
a recent review, see Ref. [19]). In a manner analogous to the proton PDFs, the nPDFs are 
needed in order to predict observables in proton-nucleus (pA) and nucleus-nucleus (AA) 
collisions. However, as compared to the proton case, the nuclear parton distribution func- 
tions are far less well constrained. Data that can be used in a global analysis are available 
for fewer hard processes and also cover a smaller kinematic range. In particular, the nuclear 
gluon distribution is only very weakly constrained, leading to a significant uncertainty in 
the theoretical predictions of hard processes in AA collisions. 
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For this reason it is crucial to use a variety of hard processes in pA cohisions, both 
at RHIC and at LHC, in order to better constrain nuclear parton densities. The inclu- 
sive production of jets, lepton pairs or vector bosons are natural candidates since they 
are already used in global analyses of proton PDFs^. In addition, other processes which 
could constrain the gluon nPDF have been discussed in the literature and have yet to be 
employed. For instance, the production of isolated direct photons [21] as well as inclu- 
sive hadrons [22] at RHIC and LHC can provide useful constraints on the nuclear gluon 
distribution^, even though in the latter channel the fragmentation process complicates its 
extraction. Another natural candidate for measuring the gluon nPDF is heavy-quark [23] 
or heavy-quarkonium [24] production. Quarkonium production is however still not fully 
under control theoretically (see e.g. [25] for a review), hence it is not obvious whether a 
meaningful extraction of the nuclear gluon PDF will eventually be possible in this channel, 
yet indirect constraints might be obtained [26]. 

In this paper, we investigate the production of a direct photon in association with a 
heavy-quark jet in pA collisions in order to constrain parton densities in nuclei^. As we 
will show, this process is dominated by the heavy-quark-gluon (Qg) initial state at both 
RHIC and the LHC making the nuclear production ratio in pA over pp collisions, 

„7Q ^ ajpA^'y QX) 

Aaipp^jQX)' ^ ' 

a useful observable in order to determine the gluon and heavy-quark nPDFs in complemen- 
tary x-ranges from RHIC to LHC. One of the advantages of such a ratio is that many of the 
experimental and theoretical uncertainties cancel. Nevertheless, for a solid interpretation of 
the ratios it is also necessary to compare the theory directly with the (differential) measured 
cross-sections. For this reason we present cross-sections and pT-distributions computed at 
next-to-leading order (NLO) of QCD using acceptance and isolation cuts appropriate for 
the PHENIX and ALICE experiments at RHIC and LHC, respectively. Using the available 
luminosity values we also provide simple estimates for the expected event numbers. 

The paper is organized as follows. In section 2 we briefly describe the NLO calculation 
used in the present paper (more details can be found in [28, 29]). In section 3, we discuss 
the different nPDF sets used in our analysis, focusing especially on the gluon and the heavy- 
quark sectors. In sections 4 and 5, results in pA collisions at RHIC and LHC, respectively, 
are presented. In each case, we start with a discussion of the acceptance and isolation cuts, 
then turn to the (differential) cross-sections and event numbers, followed by a discussion 
of the nuclear production ratios. Finally, we summarize our main results in section 6. 

^Recently, a paper by Paukkunen and Salgado [20] discussed that weak boson production at the LHC 
might be useful in order to constrain nPDFs. 

^Note that the EPS08/EPS09 [14, 15] global analyses include single-inclusive iv" data from the PHENIX 
experiment at RHIC. 

recent paper by Betemps and Machado [27] has performed a calculation for 7 + c production using 
however the target rest frame formalism. 
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2 Direct photon production in association with a heavy-quark jet 

Single direct photons have long been considered an excellent probe of the structure of the 
proton due to their point-like electromagnetic coupling to quarks and due to the fact that 
they escape confinement [30, 31]. Their study can naturally be extended to high-energy 
nuclear collisions where one can use direct photons to investigate the structure of nuclei as 
well [21]. 

However, it might also be relevant to study more exclusive final states, such as the 
double inclusive production of a direct photon in association with a heavy-quark (charm, 
bottom) jet ^ in order to get additional constraints on parton distribution functions. The 
lower counting rates expected for this observable are compensated by various advantages: 

• As shown below, the cross-section for direct photon plus heavy-quark production in 
pp and pA collisions is largely dominated by the gluon-heavy-quark {gQ) channel. 
This offers in principle a direct access to the gluon and heavy-quark distributions in 
a proton and in nuclei; 

• A two-particle final-state allows for the independent determination of the parton 
momentum fractions xi (projectile) and X2 (target), using leading order kinematics 
and in the absence of fragmentation processes; 

• Since the valence up quark distribution (to which single photons mostly couple) is 
smaller in neutrons - and therefore in nuclei - as compared to that in a proton, the 
nuclear production ratio ii^^ of single photon production at large Xj. = 2pj,/y/s is 
different than 1 independently of any nPDF effects [32]. In the 7 + Q production 
channel the photon couples mostly to the heavy-quark, which, by isospin symmetry, 
has the same distribution in a proton or neutron, i.e. = Q", leading to a nuclear 
production ratio i?^^ free of any "isospin" effects and thus properly normalized to 1 
in the absence of nPDF corrections. 

At leading-order accuracy, 0{aas), at the hard-scattering level the production of a 
direct photon with a heavy-quark jet only arises from the gQ — )■ Compton scattering 
process, making this observable highly sensitive to both the gluon and heavy-quark PDFs. 
This is at variance with the single photon channel for which the Compton scattering (gq — )• 
75) as well as annihilation process {qq — )■ ^g) channels compete^. At NLO the number 
of contributing subprocesses increases to seven, listed in Table 1. As can be seen, all 
subprocesses apart from qq — )• 7QQ are g and/or Q initiated. Which of these subprocesses 
dominate is highly dependent on the collider type {pp vs. pp/pA) and the collider center- 
of-mass energy. For example, g and Q initiated subprocesses will be more dominant at pp 
and pA colliders, whereas at the Tevatron [pp collisions) the qq — t- 7QQ dominates at high 
Pj. because of the valence-valence qq scattering in these collisions. 

*We do not distinguish between heavy-quark and heavy-anti-quark in the final state, that is the sum of 
7 + Q and 7 -I- Q is considered; therefore here by qQ — ^ qQ^ we choose to denote the sum of qQ — > qQ'y, 
qQ — ^ qQ'y, qQ qQj and qQ qQ-y. 

""In pp and pA coUisions, however, the Compton subprocess largely dominates the annihilation process 
from the dominance of the gluon distribution over that of sea-quarks. 
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99 iQQ 9Q 79Q 

Qq —?■ ^qQ Qq —?■ ^qQ 
QQ -iQQ QQ ^ ^QQ 
qq jQQ 

Table 1. List of all 2 — 3 NLO hard-scattering subprocesses. 



When one considers higher order subprocesses, such as qQ — )• qQj, the produced 
photon may be emitted collinearly with the final state q giving rise to a collinear singularity. 
This singular contribution is absorbed in fragmentation functions (FFs) D.y/g[z, fi'^), which 
satisfy a set of inhomogeneous DGLAP equations, the solutions of which are of order 
O (a/as). As a consequence, another class of contributions of order O (aag) consists of 

2 — )■ 2 QCD subprocesses with at least one heavy-quark in the final state and another 
parton fragmenting into a collinear photon. These so-called fragmentation contributions 
need to be taken into account at each order in the perturbative expansion. As in the 
LO direct channel, we also include the O (aa^) fragmentation contributions, which are 
needed for a complete NLO calculation. It should however be mentioned that isolation 
requirements ~ used experimentally in order to minimize background coming from hadron 
decays - greatly decrease these fragmentation contributions. 

The present calculations have been carried out using the strong coupling constant 
corresponding to the chosen PDF set: a^^'^(M^) = 0.118 in next-to-leading order for both 
nCTEQ and EPS09, and a^^'^(M^) = 0.1165 for HKN. The renormalization, factorization 
and fragmentation scales have been set to hr = /if = fJ-f = PT-y and we have used rUc = 1.3 
GeV and rub = 4.5 GeV for the charm and bottom quark masses. We utilize the photon 
fragmentation functions of L. Bourhis, M. Fontannaz and J. P. Quillet [33]. For further 
details on the theoretical calculations, the reader may refer to [28, 29]. 

3 Nuclear Parton Distribution Functions 

In order to obtain results in hadronic collisions, the partonic cross-sections have to be 
convoluted with PDFs for protons and nuclei. For the latter we show results using the 
most recent nCTEQ [9, 10], EPS09 [15], and HKN07 [18] nuclear PDF sets^. Each set 
of nuclear PDFs is connected to a set of proton PDFs to which it reduces in the limit 
j4 — >■ 1 where A is the atomic mass number of the nucleus^. Therefore we use the various 
nPDFs together with their corresponding proton PDFs in the calculations. Since our goal 
is to probe gluon and heavy-quark nPDFs, let us now discuss these specific distributions 
in greater detail. 

^Note that the nDS04 PDFs [12] are not considered here since these are obtained in a 3-fixed flavor 
number scheme (no charm PDF) whereas our calculation is in a variable flavor number scheme. 

■^More precisely, EPS09 is linked to the CTEQ6.1M proton PDFs [34], HKN07 to the MRST98 [35] set, 
and the nCTEQ PDFs to the reference PDFs described in Ref. [36] which are very similar to the CTEQ6.1M 
distribution functions [34]. This reference set excludes most of the nuclear data used in the PDF global flt, 
and therefore is not biased by any nuclear corrections. 
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Figure 1. Nuclear modifications R-^ = gP^^{x,Q)/gP{x,Q). Left: for gold at Q = 15 GcV. Right: 
for lead at Q = 50 GeV. Shown are results for nCTEQ decutS (solid, black line), EPS09 (dashed, 
blue line) + error band, HKN07 (dash-dotted, red line) + error band. The boxes exemplify the 
x-regions probed at RHIC (Vs^ — 200 GcV) and the LHC (y's^ — 8.8 TeV), respectively. 



3.1 Gluon sector 

As already mentioned, the nuclear gluon distribution is only very weakly constrained in 
the a;-range 0.02 < x < 0.2 from the Q^-dependence of structure function ratios in deep- 
inelastic scattering (DIS) [37], F^''{x,Q^)/F^{x,Q^), measured by the NMC collabora- 
tion [38] 8. 

In order to compare the various nPDF sets, we plot in Fig. 1 the gluon distribution 
ratio R^{x,Q) = g'P^^{x,Q)/gP{x,Q) as a function of x for a gold nucleus at Q = 15 GeV 
(left) and for a lead nucleus at Q = 50 GeV (right). The chosen hard scales Q = 15,50 
GeV are typical for prompt photon production at RHIC and the LHC, respectively, and 
the boxes highlight the x-regions probed by these colliders. 

As can be seen, the nuclear gluon distribution is very poorly constrained^, especially in 
the regions x < 0.02 and x > 0.1. The uncertainty bands of the HKN07 and EPS09 gluon 
distributions do not overlap for a wide range of momentum fractions with x > 0.02. Also 
the rather narrow and overlapping bands at small x < 0.02 do not reflect any constraints by 
data, but instead are theoretical assumptions imposed on the small-x behavior of the gluon 
distributions. The nCTEQ gluon has again quite a different x-shape which is considerably 
larger (smaller) in the x-region probed by RHIC (the LHC) as compared to HKN07 and 
EPS09. 

At present, the nCTEQ nPDFs do not come with an error band. In order to assess 
the uncertainty of the nuclear gluon PDF we have performed a series of global fits to iA 
DIS and Drell-Yan data in the same framework as described in Ref. [10]. However, each 
time we have varied assumptions on the functional form of the gluon distribution^'^. More 

*As discussed earlier, EPS08/EPS09 also include inclusive vr'^ data from the PHENIX experiment at 
RHIC, with a strong weight in order to better determine the nuclear gluon distribution. 

^Note also that at lower scales the uncertainties of the nPDFs are even more pronounced. 
^°The corresponding sets of nPDFs are available upon request from the authors. 



- 5 - 



Name 


(initial) fit parameter 


Cl,l 


Cl,2 


decutS 


free 


-U.Zy 


-u.uy 


decutSgl 


nxed 


U.2 


50.0 


decut3g2 


nxed 


-U.l 


-0.15 


decut jgj 


fixed 


0.2 


-0.15 


decut3g4 


free 


0.2 


-0.15 


decut3g5 


fixed 


0.2 


-0.25 


decutSg? 


fixed 


0.2 


-0.23 


decutSgS 


fixed 


0.35 


-0.15 


decut3g9 


fixed - free proton 


0.0 





Table 2. Start values for the parameter ci = ci o + ci_i(l — A^'^'^'^) governing the small x behavior 
of the gluon distribution at the initial scale Qo = 1-3 GeV. The parameter ci,o corresponds to the 
gluon in the proton and has been kept fixed. With one exception, decut3g4, the parameters ci^i 
and ci^2 have been kept fixed as well. For further details on the functional form the reader may 
refer to Rcf. [10]. 




Figure 2. Left: nPDF ratio R^'' at Qo = 1.3 GeV predicted within the different nCTEQ sets 
~ fits from top to bottom: decut3g9, decut3g5, decut3g7, decut3g8, decut3g3, decut3g4, 
decut3g2, decut3gl, decut3g. Right: nCTEQ gluon nPDFs for different A (1, 2, 4, 9, 12, 27, 
56, 108, 207) vs x at Qo = 1-3 GeV - from left to right, and top to bottom: decut3g, decut3gl, 
decut3g2, decut3g3, decut3g4, decut3g5, decut3g7, decut3g8, decut3g9. 

precisely, the coefficient ci = ci^o + ci,i(l — A^^'^''^) influencing the small x behavior of the 
gluon distribution, see Eq. (1) in [10], has been varied as summarized in Table 2. Each of 
these fits is equally acceptable with an excellent x^/dof in the range of x^/dof= 0.88-0.9. 

In order to give an idea about the gluon nPDF uncertainty, we plot in Fig. 2 (left) 
a collection of ratios Rg^ for a lead nucleus as a function of the momentum fraction x at 
the initial scale Qq = 1.3 GeV, while in Fig. 2 (right) the actual gluon nPDFs are plotted 
versus x for a range of A values. Results are shown for several of the fits of the decutSg 
series. The ensemble of these curves together with the HKN07 and EPS09 uncertainty 
bands provides a much more realistic estimate of the uncertainty of the nuclear gluon 
distribution which is clearly underestimated by just one individual error band. This is due 
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Figure 3. Nuclear modifications to deutcron, R'^g = gP/'^{x,Q)/gP{x,Q) at Q 15 GeV, nCTEQ 
(solid black line), EPS09 (dashed blue line), HKN (dash-dotted red line) + error band 

to the fact that for a specific fit, assumptions on the functional form of the nPDFs have 
been made so that the error bands based on the Hessian matrix for a given minimum only 
reflect the uncertainty relative to this set of assumptions. 

In order to explore the ahowed range of nCTEQ predictions for the nuclear produc- 
tion ratios to be discussed in Sec. 4 and 5 we choose the three sets decutS (solid black 
line), decut3g9 (dotted red line), and decutSgS (dash-dotted green line). The original 
fit decutS [10] exhibits a very strong shadowing at small x; conversely, the decut3g9 fit 
closely follows the distribution of the gluon in a (free) proton and the decutSgS gluon lies 
between the two extremes. In most cases, however, we focus on the original fit decutS to 
which we refer by default as nCTEQ, if the fit name is not specified. Together, with the 
HKN07 and EPS09 predictions this will cover to a good degree the range of possibilities 
for the nuclear production ratios. 

At RHIC, the incoming projectile is not a proton but a deuteron nucleus {A = 2), whose 
PDFs may be different from that of a proton. In Fig. 3 the expected nuclear modifications 
of the deuteron nucleus are shown. The EPS09 nPDFs do not include nuclear corrections 
to the deuteron PDFs, while the HKN and nCTEQ sets do. Those corrections are not 
large, at most 5%, with nCTEQ having them more pronounced. 

3.2 Heavy-quark sector 

Let us now turn to the heavy-quark distribution. In the standard approach used in almost 
all global analyses of PDFs, the heavy-quark distributions are generated radiatively, ac- 
cording to DGLAP evolution equations [39-41], starting with a perturbatively calculable 
boundary condition [42, 43] at a scale of the order of the heavy-quark mass. In other 
words, there are no free fit parameters associated to the heavy-quark distribution and it is 
entirely related to the gluon distribution function at the scale of the boundary condition. 
As a consequence, the nuclear modifications to the radiatively generated heavy-quark PDF 
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Figure 4. Left: nPDF ratios = gP^P''{x,Q)/gP{x,Q) (top) , = cP^^''{x,Q)/cP{x,Q) 

(middle), R^^ = ul^^\x,Q)/uP{x,Q) (bottom) at Q = 50 GeV witliin nCTEQ (solid black line), 
EPS09 (dashed blue line), and HKN07 (dash-dotted red line). The shaded regions correspond to 
the x-values probed at RHIC {x - lO^i) and the LHC (a; ~ lO^^). Right: double ratios -RfV-R^ 
and R^^'/R^l using the same nPDF sets. 



are very similar to those of the gluon distribution and quite different from the nuclear 
corrections in the valence-quark sector. This feature is illustrated in Fig. 4 (left) where 
we show the nuclear modifications for the gluon (upper panel), charm (middle panel) and 
the valence up-quark (bottom panel), in a lead nucleus, for three different sets of nuclear 
PDFs at the scale Q = 50 GeV as in Fig. 1 (right). The shaded regions in Fig. 4 (left) 
correspond to the typical x- values probed at RHIC (x ~ 10~^) and the LHC (x ~ 10~^). 
The close similarity between the charm and the gluon nPDFs can be better seen in Fig. 4 
(right) where the double ratios, R^^ /Rg^ and R^^/R^^ {u^, = u — u being the valence 
distribution), are plotted. Remarkably, the nuclear effects in the gluon and the charm 
PDFs are different by at most 20% at large x {Rc/Rg ^ 1-2), whereas the difference can 
be as large as 80% {Rc/Ruy — 1-8) when comparing the valence up-quark and the charm 
nPDF ratios. Therefore, in the standard approach, the LO direct contribution (gQ — t- 7Q) 
only depends on the gluon distribution, either directly or via the dynamically generated 
heavy-quark distribution, making this process an ideal probe of the poorly known gluon 
nPDF. 

Conversely, light-cone models predict a nonperturbative (intrinsic) heavy-quark com- 
ponent in the proton wave-function [44, 45] (see [46] for an overview of different models). 
Recently, there have been studies investigating a possible intrinsic charm (IC) content in 
the context of a global analysis of proton PDFs [1, 47]. In the nuclear case, there are 
no global PDF studies of IC (or IB) available. This is again mainly due to the lack of 
nuclear data sensitive to the heavy-quark components in nuclei. For this reason, we only 
consider the standard radiative charm approach in the present paper. Measurements of 

"Note that in Melhn moment iV-space the relation c''^'^{N,Q)/c''{N,Q) ~ g^'^[N ,Q) / {N ,Q) holds 
approximately for Q ~ Qo- 
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7 + Q production in pA collisions at backward (forward) rapidities are sensitive to the 
BHPS-IC in nuclei (the proton) complicating the analysis. A similar statement is true for 
RHIC, where due to the lower center-of-mass energy the results depend on the amount of 
intrinsic charm. Therefore, once the nuclear gluon distribution has been better determined 
from other processes these cases may be useful in the future to constrain the nuclear IC. 

4 Phenomenology at RHIC 

In this section we present the theoretical predictions for the associated production of a 
photon and a heavy-quark jet in d-Au collisions at RHIC at ^/s^ = 200 GeV. 

4.1 Cuts 

The experimental cuts used for the theoretical predictions are listed in Table 3. The photon 
rapidity and isolation requirements are appropriate for the PHENIX detector [48]. When 
Pt,7 = Pt,q the NLO cross-section is known to become infrared sensitive^^. Therefore, 
in order to acquire an infrared safe cross-section, the minimum transverse momentum of 
the photon is kept slightly above that of the heavy-quark [49, 50] which ensures a proper 
cancellation between real and virtual contributions. Also note that the p^*" cuts in the 
J + b channel (p^q = 14 GeV and = 17 GeV) were taken to be higher than those in 
7 -|- c events in order to keep terms of O [rnq/pT) small. 





Pt 


Rapidity 




Isolation Cuts 


Photon (-|-c) 


Pt,^ 


= 7 


GeV 


\y^\ < 0.35 


0° < (J 


b < 180° 


R = 0.5, 


e < O.IE^ 


Photon {+h) 


Pt,"i 


= 17 


GeV 


\y^\ < 0.35 


0° < (J 


b < 180° 


R = 0.5, 


e < O.IE^ 


Charm Jet 


min 

Pt,q 


= 5 


GeV 


lygl < 0.8 










Bottom Jet 


Pt,q 


= 14 


GeV 


IvqI < 0.8 











Table 3. Experimental cuts used for the theoretical predictions at RHIC. 



4.2 Spectra and expected rates 

The spectra are shown for -y + c production in Fig. 5 (left) and 7 -I- 6 production 
in Fig. 5 (right) where the band represents the scale uncertainty obtained by varying 
the renormalization, factorization and fragmentation scales by a factor of two around the 
central scale choice, i.e., /ir = hf = fJ-f = CpT'j with ^ = 1/2, 2. 

The total integrated cross-section for 7 -|- c events is cri^+c = 37036 pb. Using the 
projected weekly luminosity for dAu collisions at RHIC-II, C^^^^ = 62 nb^^ [51], and 
assuming 12 weeks of ion runs per year, the yearly luminosity is £2'^"''' = 744 nb~^. Thus, 
an estimate of the number of events expected in one year is Nj^j^^ = ly^"''^ x cri^^" ~ 2.8 x 10"^ 
in d-Au collisions, without taking into account effects of the experimental acceptances and 
efficiencies. At pT^ ~ 20 GeV (da/dpT^ ~ 45 pb/GeV), the number of events would still 

^^This back-to-back kinematics matches the LO case and constrains the transverse momentum of the 
third particle to be zero. 
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Figure 6. Subprocess contributions to the differential cross-section at RHIC: NLO (solid black 
line), LO+gg QQj (dashed blue line), gQ — )• gQ-f (dash-dotted purple line), — > QQ'j (dotted 
red line), q{q)Q{Q) q{q)Q{Q)j,QQ QQl (dash-dot-dottcd magenta Hne). 



be large, O (lO^) per GeV-bin. This indicates that the number of 7-l-c events in a year 
produced at RHIC-II will be substantial. The rates expected in the 7 -|- 6 channel at RHIC 
are naturally much more modest. Using the total integrated cross-section o"^^^ = 32 pb, 
the number of events to be expected in a year is N^^^ = 24. Therefore we shall mostly 
focus the discussion on the 7-l-c channel in the following. 

In Fig. 6 the individual subprocess contributions to the 7-l-c NLO production cross- 
section are presented. As can be seen, the dominant subprocesses are the LO Compton 
scattering gQ — > ^Q, as well as the higher-order gQ — t- 'jgQ and gg — )■ 'jQQ channels. 
Thus almost all the PDF dependence in the NLO 7-l-c cross-section comes from the gluon 
and heavy-quark PDFs and not from the light-quark PDFs. The relative increase of the 
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Figure 7. Left: nuclear produetion ratio of the 7 + c cross-section at RHIC using nCTEQ (solid 
black line), nCTEQ without nuclear corrections in the deuteron (dotted magenta line), EPS09 
(dashed blue line) + error band, HKN (dash-dotted red line) -I- error band. Right: nuclear mod- 
ification of the gluon in gold, Rf^{x,Q = xy/S/2 - pt), for the x- region probed at RHIC. This 
figure corresponds to the enlargement of the box region in the left panel of Fig. 1 . 
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Figure 8. Nuclear production ratio of the 7 -f c cross-section at RHIC using the three nCTEQ 
fits discussed in Sec. 3. Also shown is the scale dependence of 



contributions by the annihilation subprocess, qq — )• QQ'J, and the light quark-heavy quark 
subprocess qQ — )• qQj at higher x {pT-y ~ 15 GeV) is due to the slower decrease of the 
valence quark PDF at high x as compared to the rest of the PDFs. 



4.3 Nuclear production ratios 

Let us now discuss the nuclear modifications of 7 -|- c production in d-Au collisions. The 
nuclear production ratio, 



o7+c 
^dAu 



1 da / dpT-yidkn ^ 'y + c + X) 
2x197 da/dpT-yipp J + c + X) ' 



(4.1) 
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is plotted in Fig. 7 (left) as a function of px-y using the three nPDF sets discussed in 
section 3, namely nCTEQ (solid black line), EPS09 (dashed blue line + error band) and 
HKN (dash-dotted red line + error band). There is some overlap between -R^kn -^epsoq 
at not too large < 15 GeV, whereas the difference between -RnCTEQ hand 
and -Rhkn -^epso9 other hand is larger for all transverse momenta. The -RnCTEQ 

ratio is further increased by the anti-shadowing corrections in the deuteron projectile, as 
can be seen in Fig. 7 (left) where the nCTEQ predictions are performed with (solid line) 
and without (dashed) corrections in the deuteron (see also Fig. 3). Due to the rather low 
center-of-mass energy (as compared to the Tevatron/LHC) the collisions at central rapidity 
at RHIC probe relatively high values of momentum fractions carried by the partons in the 
nuclear target, = 0{2pr^/y/s) = 0(10"^). In Fig. 7 (right) we show the nuclear 
modifications of the gluon distribution in a gold nucleus, R^^{x,Q = x^/S/2 ~ px), for 
the typical x-region probed at RHIC. Note that, this figure corresponds to the enlargement 
of the box-region in the left panel of Fig. 1. As can be seen the nuclear production ratios of 
7 + c events shown in Fig. 7 (left) closely correspond to the different nuclear modifications 
of the gluon distribution depicted on the right side of Fig. 7. Clearly, measurements of 
this process with appropriately small error bars will be able to distinguish between these 
three different nuclear corrections to the cross-section and therefore be able to constrain 
the gluon nuclear PDF. 

In Fig. 8 we present the dependence of the nuclear modifications on the three nCTEQ 
fits (decutS, decut3g9, decutSgS) discussed in Sec. 3. It is clear that these different fits 
cover quite a spread of nuclear modifications, ranging from ones which are quite pronounced 
(decutS) to almost none (decut3g9). We stress again, that neither of these predictions is 
preferred over the other since the nuclear gluon distribution is so poorly known. Finally, 
we also show in Fig. 8 the scale uncertainty which is entirely negligible compared to the 
PDF uncertainty. 

In the next section we present the phenomenology of 7 + Q production at the LHC 
where smaller values of x^ are probed due to the higher center-of-mass energy. 



5 Phenomenology at LHC 

In this section, calculations are carried out for p-Ph collisions at the LHC nominal energy, 
^/s^ = 8.8 TeV, different from the pp collision energy {^/s = 14 TeV). 

5.1 Cuts 

The cuts used in the present calculation are shown in Table 4 and are appropriate for 
the ALICE detector^^ [52-56]. Note that the rapidity shown in Table 4 is given in the 
laboratory frame, which in pA collisions is shifted by Ay = —0.47 with respect to the center- 
of-mass frame [57, 58]. In ALICE, photons can be identified in the EMCal electromagnetic 
calorimeter, or in the PHOS spectrometer with a somewhat more limited acceptance. 

^^We have verified that similar results and conclusions are obtained when using either ATLAS or CMS 
acceptances at central rapidities. 
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PT 


Rapidity 




Isolation 


Cuts 


Photon (PHOS) 




= 20 


GeV 


\y^\ < 0.12 


220° <<p< 320° 


R = 0.2, 


„th 

Pt 


= 2 GeV 


Photon (EMCal) 


^min 


= 20 


GeV 


\y^\ < 0.7 


80° <(j)< 180° 


R = 0.2, 


pf 


= 2 GeV 


Heavy Jet 


Pt,q 


= 15 


GeV 


IvqI < 0.7 











Table 4. Experimental cuts for the ALICE detector. 



5.2 Spectra and expected rates 

The differential NLO cross-section is plotted as a function of the photon transverse mo- 
mentum in the 7 + c (7 + &) channel in Fig. 9 left (right) for both PHOS (lower band) and 
EMCal (upper band); the dotted curves indicate the theoretical scale uncertainty. 

In order to estimate the number of events produced, we use the instantaneous lumi- 
nosity = 10"'' pb~^s~^ [58] which corresponds to a yearly integrated luminosity of 
£year _ iQ-1 pb^^ assuming one month (At = lO^s) of running in the heavy-ion mode at 
the LHC. In Table 5 the total integrated cross-section for 'j + Q for both PHOS and EMCal 
along with the respective anticipated number of events (before experimental efficiencies), 
^j+Q ~ ^ are given. As expected the 7+6 and 7+c cross-sections at EMCal are 
increased substantially by the larger acceptance of that detector. The number of expected 
7 + 6 events is large, at variance with what is expected at RHIC (see section 4.2). 





pPb 




7 -he PHOS 


22700 pb 


2270 


7 + 6 PHOS 


3300 pb 


330 


7 + c EMCal 


119000 pb 


11900 


7 + 6 EMCal 


22700 pb 


2270 



Table 5. Total integrated cross-section and number of events per year for 7 + Q production in 
p-Pb collisions at the LHC for PHOS and EMCal acceptances. 



The individual subprocess contributions to the cross-section are depicted in Fig. 10. 
As one can see the Compton {gQ — )• jQ) as well as the gQ — )• jgQ and gg — >• jQQ are 
the dominant subprocesses, demonstrating the sensitivity of this process to the gluon and 
charm nPDFs. Here the contribution by the annihilation subprocess proves much smaller 
than at RHIC. This is caused by the less pronounced difference in the light anti-quark 
and heavy quark PDFs at small x as compared to large x. So that now qq — t- ^QQ can 
no longer compete with the light quark-heavy quark (antiquark) piece of the cross-section 
{qQ qQ-f, qQ qQ-/, qQ qQ-f, qQ qQ-f). 

5.3 Nuclear production ratios 

The nuclear production ratio R'^p^ = ^ '^di'/ZTXp^i+c+x) shown in Fig. 11 (left) us- 
ing the nCTEQ decutS (solid black line), nCTEQ decut3g3 (dotted black line), nCTEQ 
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Figure 9. NLO differential cross-section for 7 + c (left) and 7 + 6 (right) production in p-Pb 
collisions at a center-of-mass energy of y/s^ = 8.8 TeV in PHOS (lower band) and EMCal (upper 
band) acceptances. 




50 100 150 



Px/GeV) 

Figure 10. Subprocess contributions to the differential cross-section shown in Fig. 9 (left), NLO 
(soHd black line), LO+gg — )- QQj (dashed blue line), gQ — >■ gQj (dash-dotted purple Hne), qq — )• 
QQj (dotted red line), qQ — >■ qQ^; QQ — >■ QQ7 (dash-dot-dotted magenta line). 

decut3g9 (dash-dot-dashed black line), EPS09 (dashed blue line), and HKN07 (dash- 
dotted red line) nuclear PDFs. For the latter two cases the bands represent the nPDF 
uncertainties calculated as described in section 3. Remarkably, there is almost no overlap 
between the EPS09 and the HKN predictions, therefore an appropriate measurement of 
this process will be able to distinguish between the two nPDF sets. The nCTEQ nuclear 
modification, using the decutS fit, is considerably different from the two other sets at lower 
values of Pj,. On the other hand, for the decut3g9 set, the nuclear modification factor is 
close to unity in the x-range as shown in Fig. 11 (right), giving rise to the nuclear pro- 
duction ratio for this nCTEQ set which lies inside the EPS09 uncertainty band. Fig. 11 
(left). We stress again that both, decutS and decut3g9, are perfectly acceptable fits to 
the iA DIS-I-DY data with different assumptions on the small x behavior. We further 
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PT.(GeV) 



Figure 11. Left: nuclear production ratio of 7 + c cross-section at LHC within ALICE PHOS 
acceptances, using nCTEQ decutS (solid black line), nCTEQ decutSgS (dotted black line), nCTEQ 
decut3g9 (dash-dot-dashcd black line), EPS09 (dashed blue line) + error band, HKN07 (dash-dotted 
red line) + error band. Right: R^''{x, Q = x\/S /2 ^ pr) ratio as a function of x, in the x region 
probed at the LHC. This figure corresponds to the enlargement of the box region in the right panel 
of Fig. 1. 

show the ratio for decut3g3, as a representative lying between the two extremes. Inspect- 
ing Fig. 2, it is clear that the rest of the predictions from the decutSg series would fill 
the gap between the decutS curve and the decut3g9 curve. Taken together, this gives 
a more realistic impression of the true PDF uncertainty of the nuclear production ratio. 
Therefore, measurements in this region will provide useful constraints on the nuclear gluon 
distribution. 

Some further comments are in order: (i) In this paper we have demonstrated that 
the ratio of the j + c cross-section in pA over pp collisions at central rapidities will be 
very useful to constrain the nuclear gluon distribution. At forward rapidities, even smaller 
values could be probed in the nuclear targets where the uncertainties are largest. At 
backward rapidities, large is probed, hence the cross-section in this rapidity region will be 
sensitive to any existent intrinsic charm contribution in the nucleus. Such a measurement 
could be performed with the CMS and ATLAS detectors which cover a wider range in 
rapidity. We postpone such a study to a future publication, since currently there are 
no available IC nuclear PDFs; (ii) At the LHC ^ + b events will also be produced with 
sufficient statistics. Experimentally this channel might be preferable due to the much 
better 6-tagging efficiencies. Furthermore, uncertainties related to possible intrinsic charm 
contributions should be much reduced in the bottom case. However, as for j+c production, 
the nuclear production ratios follow closely the gluon ratio and, therefore, we do not show 
a separate figure here. 
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6 Conclusions 



We have performed a detailed phenomenological study of direct photon production in 
association with a heavy-quark jet in pA cohisions at RHIC and at the LHC, at next-to- 
leading order in QCD. The dominant contribution to this process is given by the gQ — )• 
7Q[+5] subprocess. This offers a sensitive mechanism to constrain the heavy-quark and 
gluon distributions in nuclei, whose precise knowledge is necessary in order to predict the 
rates of hard processes in heavy-ion collisions where quark-gluon plasma is expected to be 
formed. 

We have performed the calculation of 7 + Q production spectra at RHIC and at the 
LHC within the acceptances of various detectors (PHENIX and ALICE-PHOS/ALICE- 
EMCal) and have presented the corresponding counting rates. At the LHC the j + c and 
7 -|- 6 production rate is important, while at RHIC only 7 -|- c events will be copiously 
produced. 

Our results for RHIC (see Fig. 7) exhibit a strong sensitivity to the nuclear gluon 
distribution permitting to constrain it at a; ~ 0.1-0.2. Similarly to RHIC the ratio at the 
LHC (see Fig. 11) is very sensitive to the gluon distribution probing a smaller x ~ 10~^, 
i.e. in a complementary range to RHIC. These results have been obtained in the "standard 
approach" of radiatively generated charm distribution. A future study will focus on the 
possibility to constrain the intrinsic charm contribution to the nucleus as well as the proton. 
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